JOURNAL OF XXX, VOL. X, NO. X, DATE GOES HERE 



1 
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Abstract — Link adaptation is a crucial part of many modern 
communications systems, allowing the system to adapt the trans- 
mission and reception strategies to changes in channel conditions. 
One of the fundamental components of the link adaptation 
mechanism is signal to noise ratio (SNR) estimation, measuring 
the instantaneous (mostly post processing) SNR at the receiver. 
That is, the SNR at the decoder input, which is an important 
metric for the prediction of decoder performance. In linearly 
decoded MIMO, which is the common MIMO decoding strategy, 
the post processing SNR is well defined. However, this is not 
the case in optimal maximum likelihood (ML) decoding applied 
to spatial multiplexing (SM). This gap is interesting since ML 
decoded SM is gaining ever growing interest in recent research 
and practice due to the rapid increase in computation power, 
and available near optimal low complexity schemes. In this paper 
we close the gap and provide SNR estimation schemes for ML 
decoded SM, which are based on various approximations of 
the "per stream" error probability. The proposed methods are 
applicable for both horizonal and vertical decoding. Moreover, 
we propose a very low complexity implementation for the SNR 
estimation mechanism employing the ML decoder itself with 
negligible overhead. 

Index Terms — SNR Estimation, CINR, MIMO, Spatial Multi- 
plexing, Maximum Likelihood Decoding. 



I. Introduction 

ESTIMATING the signal to noise ratio (SNR) is one of the 
important tasks in communications systems. The measure 
of SNR indicates the quality of the channel and enables the use 
of link adaptation to improve the spectral efficiency. The main 
idea of link adaptation is to use the transmission parameters 
that yield the highest possible bit rate. The most common 
parameters to be adapted are the modulation and coding 
scheme. In addition some other parameters may be adjusted 
for the benefit of the systems such as transmit power levels, 
bandwidth usage or MIMO mode (when a MIMO scheme is 
applied). 

Extensive work has been performed on this area and several 
approaches were introduced. In recent years, a number of new 
link adaptation schemes have been proposed for different types 
of wireless networks. An example for that is a Receiver-Based 
Auto-Rate (RBAR) protocol based on the RTS/CTS (Request- 
To-Send/Clear- To-Send) mechanism (TJ. The basic idea of 
RBAR can be summarized as follows. First, the receiver 
estimates the wireless channel quality using a sample of the 
instantaneously-received signal strength at the end of the RTS 
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reception. The receiver selects the appropriate transmission 
rate based on this estimate, and feeds back to the transmitter 
using the CTS. Then, the transmitter responds to the receipt 
of the CTS by transmitting the data packet at the rate chosen 
by the receiver. 

In the General Packet Radio Service (GPRS) development 
of GSM two link adaptation schemes were proposed [2|. One 
is based on the estimate of the Carrier to Interference ratio 
(C/I), and the other is based on the observation of the block 
error rate. 

High PErformance Radio Local Area Network type 2 
(HIPERLAN/ 2) is another wireless broadband access system 
that has been specified by European Telecommunications 
Standards Institute (ETSI) project BRAN (Broadband Radio 
Access Network) |£3j . 

Link adaptation is one of the key features of HIPERLAN/ 2 
as it has a PHY that is very similar to 802. 11a. Lin, Malmgren 
and Torsner studied the system performance of link adaptation, 
which uses the C/I as the wireless link quality measurement, 
for packet data services within HIPERLAN/2 (4)- Further- 
more, Habetha and Calvo de No presented a new algorithm 
for adaptive modulation and power control in a HIPERLAN/2 
network 0. It first assumes the maximum transmit power, and 
uses the C/I observed at the receiver to determine the proper 
PHY mode for the next frame transmission to meet the target 
packet error rate (PER). Then, it reduces the power as much 
as possible while meeting the target PER. 

A different approach uses an Euclidean distance metric to 
obtain channel quality information in terms of the average 
signal to noise ratio (SNR) (5). Then, a rate adaptation 
scheme which uses this metric to change the modulation at 
the transmitter has been described. 

Link adaptation may also be applied with a user selection 
mechanism. In user selection, we assume that multiple users 
exist and their number is larger than the number of antennas 
at the BS. This means that the BS cannot receive or transmit 
from/to all of the concurrently, so some selection mechanism 
for the formation of groups is needed. The selection process 
is to be followed by link adaptation, i.e., after the BS decides 
which group to transmit to, it should decide the modulation 
(and coding scheme in coded systems) for each user Q, (8). 

The IEEE 802.16 standard introduces 2 types of carrier to 
interference plus noise ratio (CINR) mechanisms |9|. One is 
the physical CINR (PCINR) which estimates the CINR or 
post processing CINR. The second mechanism is the effective 
CINR (ECINR) which make use of the per-tone CINR, and 
aims at proper weighting of the per-tone CINR to obtain a 
measure for the BER [ 10 1 . We emphasize that in MIMO 
transmission-reception schemes, the post processing CINR 
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(either physical or effective) is the interesting measure. Thus, 
we focus in the following sections on post processing SNR. 

II. Preliminaries 

A. SNR in MIMO Schemes 

In order to illustrate the meaning of antenna SNR (the SNR 
measured at the input to the receiver) and post processing 
SNR, we survey linearly decoded transmission and reception 
schemes. We begin with the simplest single-input single-output 
(SISO) digital communication system, having a single transmit 
(Tx) and receive (Rx) antenna. The baseband representation of 
the system is illustrated in Figure Q] The received signal y is 
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Fig. 1. SISO system block diagram 

given by 

y = hs + pv, 



(1) 



where h is the channel response which is assumed to be 
known at the receiver (perfect channel knowledge), s is the 
transmitted data (e.g QAM) and pv is an AWGN with standard 
deviation p. The instantaneous antenna SNR in this case is 
trivial and equals to the ratio between the signal's power to 
the noise's power, meaning 

r,2 ' 



SNR 



(2) 



We continue with linearly decoded MIMO systems. A multiple 
input multiple output (MIMO) system using Mr transmit 
antennas and Mr receive antennas is illustrated in Figure [2] 
The received signal y is given by 




Fig. 2. MIMO system block diagram (to be more detailed) 



y = H PHY x + pv, 



(3) 



where Hphy is an Mr x Mt channel matrix, x is the Mt x 1 
transmitted vector, v is a vector of i.i.d zero-mean complex 
Gaussian entries with unit variance and p is the noise intensity. 



In spatial multiplexing (SM) for instance, independent in- 
formation streams are transmitted through the transmit an- 
tennas. Here, the transmitted vector is x = s, where 



s = [so, s\, . . . , sm t „i] t is a vector of Mj independent 
symbols. The factor — = 



is introduced in order to maintain 



unity transmission power. The mathematical model for the 
received signal is 

V = —7==H PHY s + pn, 



(4) 



H 



where Hphy is the MIMO Mr x Mt channel matrix and 
Mr > Mt- Applying Zero-Forcing (ZF) detection gives 

s = H^y = s + pH^v = s + pGv, (5) 

where G — H* is the left pseudo-inverse of H. The post- 
processing SNR of the i— th stream in this case reads 

1 



PPSNR, = 



Si) 



1 



(6) 



p 2 (|G i>0 | 2 + ... + |G i ,M R -i| 2 )' 
Similar results are obtained for Rx diversity and Alamouti's 
space-time-coding (STC) ifTTI . In these schemes linear decod- 
ing is optimal so they may be viewed as a particular case of 

Since in general, SM is not an orthogonal transmission 
scheme, linear decoding is not optimal. The optimal ML 
decoder in this case implies exhaustive search [12|. In terms 
of post processing SNR, ML decoded SM differs from the 
linear decoding methods surveyed, in the sense that at no point 
in the reception process, there exists an expression for the 
post processing SNR. This means that another method is to 
be invoked. 

B. Known Approaches to SNR estimation in ML decoded SM 

In this section we consider some known approaches for SNR 
estimation in ML decoded SM and related issues. 

1) SNR Estimation Based on the Capacity Formula: This 
approach, that was adopted by the IEEE 802. 16e standard 
ll9l , uses capacity computation in order to estimate the SNR. 
The basic idea is to use Shannon's expression for capacity to 
evaluate the SNR, 



C 



SNR Estimate = 



1 

M 



log e det J + 



1. 



HH 



(7) 



where M is the number of independent streams. 

A natural question is the relevance of this metric to SNR. 
In order to answer this, we note that in the case of SISO, 
Rx diversity, and STC, this metric coincides with the standard 
post processing SNR. For example if we take the 1x2 MRC 
scheme where H = [ ho hi ] we obtain 



C = log 1 
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which immediately implies that the estimated SNR coincides 
with the classical ppSNR in MRC 

M 2 + M 2 



1 = 



(8) 



However, despite the fact that the ability of (0 to capture the 
SNR in the ML decoded SM case is not established, we note 
that (|7]i provides a single metric, and is thus inadequate for hor- 
izontal MIMO. Horizontal transmission indicates transmitting 
multiple separately streams over multiple antennas such that 
the number of streams is more than 1 (in contrary to vertical 
transmission that indicates transmitting a single stream over 
multiple antennas). 

2) SNR Estimation based on Error Probability Computa- 
tion: SNR and error probability are linked. Thus, an expres- 
sion for the error probability may be exploited to obtain an 
estimate for the post processing SNR (as we will demonstrate 
in the next chapter). Accordingly, we consider here the closely 
related problem of error probability calculation in ML decoded 
SM. 

The error probability in ML decoded SM does not have an 
analytic solution. One of the most prominent approaches to 
approximate and bound this probability is that of Paulraj and 
Heath lfT~3l . They obtained the following expression for the 
error probability 



III. Proposed Method for SNR Estimation in ML 
Decoded SM 

A. A Series of Approximations for the Per-Stream SNR 

We base our SNR estimation method on the asymptotic 
evaluation of the per stream error probability in ML decoded 
SM in high SNR. We emphasize that in order to obtain a 
meaningful SNR metric, the SNR estimate should satisfy (in 
QPSK) Ifl4l 

SNR, 

p(error in stream i) ss e 2 (12) 

We are therefore left with the problem of evaluating the per 
stream error probability. The conditional probability of error 
given the transmitted vector s is (throughout we condition the 
probabilities on the channel matrix H) 



Pr{error|s} = Pr (J J(S) < J(s) 

SGA(S) 



(13) 



Pr {error | if} < Q 



\d 2 . (H) 
p2 min v ' 



where d 2 ■ (H) is 
min v ; 



d 2 ■ 
"mm 



(H) = min||JT(s - c)|| 2 . 
s : ceQAM MT ,s^c 



(9) 



(10) 



where J(s) = ||Y" — Hs\\ 2 , A(s) is the set of all vectors 
corresponding to a certain type of error, such as error in 
the z — th stream which will be represented by the set Ai(s). 
Equation Qj] may be bounded using the union bound 

Pr{error|s} < £ Pr {J(3) < J(s)\a}, (14) 
seA(s) 

which may be rewritten as 

Pr {errors < £ Q 



(15) 



Since the computation of (M~) requires an exhaustive 
search, upper and lower bounds on the minimum Euclidean 
distance were introduced. 



A; 



d 2 

I tj \ min_QAM , ,2 

AM) — < cL 



d 2 

< Amaxl-n) T7 C 1 



Mt - minv / - > Mt 

These bounds depend on the maximum and minimum singular 
values of H. In case the condition number of H is high , 
these bounds may be loose. In addition, the computation of 
the singular values of H, especially for high ranked H, is a 
tough task by itself. Moreover, like the former approach, this 
method results an average estimate for the SNR over all inputs 
(streams), therefore it is not suitable for the per stream SNR 
estimation. 

The rest of the paper is organized as follows. Section II 
includes the derivation of a series of approximations for the 
SNR in ML decoded SM based on a series of approximations 
for the per stream error probability. The chapter concludes 
with a low complexity implementation of the SNR estimation 
mechanism, based on the ML decoder itself. In Section III we 
present simulation results revealing the performance (in terms 
of SNR estimation error) of the various approximations in the 
horizontal and vertical cases. In the vertical case we compare 
our result with that of standard methods. We further show that 
the QPSK based methods are valid for 16QAM and 64QAM. 
In Section IV we summarize the results and present topics for 
further research. 



seA(s) 

Focusing on events of error in the i— th stream and averaging 
w.r.t s gives 

Pr {error in i-th stream} < (16) 



Q Mt 



V2p 



seQAM t seAi(s) 

where Q~ Mt is a normalizing factor required due to the 
summation over all QAM points and transmit antennas. Using 
the upper bound on the Q-function 

Q( x )<lexp(~y (17) 

we obtain a looser but simpler bound for the probability of 
the error 

Pr {error in i-th stream} < (18) 

\\H(s-s)\\ 2 
I— V V e V 

SeQAM M T SEAi(S) 

which may be rewritten as 

Pr {error in i-th stream} < 



(19) 



1 1 



2 Q Mt 



E E - 

seQAM M T eefl»(s) 



\He\\ 2 
Ap 2 
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where e = s — s and Bi(s) is the set of vectors e corre- 
sponding to Ai(s). In order to illustrate the construction of 
the sets -Bj(s) we assume that the following symbol vector 
was transmitted 



-l+J 
V2 



(20) 



as depicted in Figure [3] We focus on the probability of error in 
the first element, i.e., the vectors e in which the first element is 
nonzero. All the possibilities for e are illustrated in Figure[3]by 
the paths from the transmitted signal to any possible erroneous 
point. These possibilities (vectors) correspond to the following 
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Fig. 3. 2 streams constellation map 

vectors e in the set Bq(s). 

B (s) = V2 



1 




1 




1 









1 




J 1 — L 


) 



1 

1-3 



-3 






-3 




-3 






■3 















1 




-1 




1 - 


- 3 




1 




1-3 ' 




' 1- 


-3 




' 1 - 


-3 






1 


~3 











1 






1 


1 




1 


~3 



(21) 



Next, we further simplify the expression, by canceling the 
dependence of the sets Bi(s) on the transmitted vector s and 
turn to unified sets Bi 



Bi = \jBi(a) 



(22) 



Using the unified sets, the error probability is approximated 

by 



\\He\\ 

Pr {error in i-th stream} w e ^P 2 

eeBi 



(23) 



Note that the approximation in d23l is based on the fact that 
elements in the sets Bi are unique, so the summation in (l23l is 
different than that in ( fT9] i. Focusing on the high SNR regime, 
(1231 may be well approximated by the max-log approximation 

as 



Pr {error in i-th stream} 



- mm — —r— 

, ees, 



This means that the sets Bi may be replaced with abbreviated 
sets Bi that omit elements that lead identical value of ||ife|| 2 , 
such that 



Pr {error in i-th stream} 



mm 

eeBi 



mm 

eeBi 



\He\\ 2 

v 

\Hef 



(25) 

T 



For instance, the vectors [ \/2 ] and [ 
lead to the same cost, so one of them may be omitted 
Moreover, noting that some pairs a, b of elements in Bi 



aa. a > 1 



(26) 



leads to the understanding that the cost of b is a 2 the cost of 
a, so b may also be omitted from the abbreviated set Bi. An 
example to such a pair is 

v/2 



V2 



and 



V2 + jV2 
V2 + jV2 

These arguments lead to the abbreviated set B 



(27) 



(28) 
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~ 1+3' 




' 1+J " 


1 




-1 




3 




—j 



An interesting outcome resulting from the abbreviation pro- 
cedure is that the sets Bi, Bj, i ^ j differ only in a single 
vector, denoted here as the first. The first vector bi in the i— th 
set, refers to the error vector in which only the i— th element 
is nonzero. This implies that the sets Bi may be written as 

Bi = {bi,C}, (30) 

where the set C is the intersection of all sets Ba, 



C = f]B l 



(31) 



(24) 



B. Low Complexity Implementation 

The structure of (|25T > implies that the SNR may be evaluated 
using the ML decoder. Such a decoder is depicted in Fig. @] 
The ML decoder uses the inputs Y and H and performs a 
search over all vectors from a predefined set A (usually A C 
QAM Mr ). The output of the decoder is mm \\Y - Hs\\ 2 . 

We notice that with setting Y = and searching over B^ 
instead of A we can evaluate the SNR introduced in d25l l. 
Since the sets Bi are significantly smaller that A and since the 
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Fig. 4. Maximum likelihood decoder 



SNR calculation process is done once for many information 
bits, the SNR estimation task embodies a small to negligible 
fraction of the ML decoder resources. 

To set ideas straight, we assume we have a 10 symbol QPSK 
allocation of 2 antenna SM. The data decoding process invokes 
the ML decoder twice for each transmitted bit. Each search 
is over ^QPSK 2 or 8 constellation points. All in all we have 
80 searches over 8 constellation points. In this allocation a 
single SNR calculation is needed, so we have the ML decoder 
invoked twice (once for each stream) to search a set of 9 
elements. This means that in this cas^e the extra overhead 

introduced by the SNR calculation is — = 2.5%. 

80 

We note that in many receivers the optimal ML decoder 
considered here is replaced with lower complexity near op- 
timal decoders (|15|,[16|). In this case, the SNR estimation 
methods may be modify to accommodate near optimal ML 
decoders. 

IV. Simulation Results 

In this chapter we present simulation studies investigating 
the performance of the series of error probability approxima- 
tions presented in the previous chapter and the corresponding 
estimated SNR. The simulation results reveal the following 
interesting virtues of the proposed SNR estimation methods. 

• The performance of the proposed methods are consistent 
with the level of approximation. This means that the 
most accurate approximation ( TTol l relying on solely on the 
union bound is superior to all other, while the simplest 
( f25l ) is inferior to all others. 

• The approximations introduce small performance degra- 
dation, so the performance of all methods is similar. 

• In case of vertical encoding (where other methods are 
applicable) the proposed methods slightly outperform the 
existing method. This serves as a good "sanity check" for 
the methods. 

• SNR should not depend on the modulation employed. 
Thus, we examine the application of a QPSK based SNR 
estimation mechanism also for higher modulations. We 
show that the QPSK based mechanism gives plausible 
results also for 16QAM. 

A. The Simulation Setup 

In order to investigate the performance of the proposed 
methods, we used the following procedure. 



1) Randomly generate 2000 i.i.d Rayleigh distributed 2x2 
matrices H. 

a) Using each matrix transmit 10 6 vectors s E 



QAM", and decode the vectors s from 
1 

y = ~^ Hs + *° n ' 



(32) 



using the ML decoder. 

b) Compute the empiric SER in each stream. 

c) Compute the empiric post processing SNR implied 
by the per stream SER according to 



SNR, = -2 log SER, 



(33) 



d) Compute the estimated SNR for QPSK at each 
stream according to the three approximations ( fTo*! ), 
d23l and (|25]) (denoted as "Union Bound", "Full- 
sum approx" and "Max-log approx" respectively in 
the graphs). 

e) Compute the empiric joint SER for both streams. 

f) Compute the empiric post processing SNR implied 
by the joint SER. 

g) Compute the joint estimated SNR for both streams 
based on the per stream estimated SNR through 
the equation 



SNR 2 [dB] , (34) 



SNR, ert [dB] = ~ (sNRi[dB] 

where SNR^ is the estimated SNR in the z — th 
stream based on ( TTSI l. (l23l and (l25l >. 
h) Compute the capacity based joint SNR, SNR cap , 
according to the capacity method (|7]i (denoted as 
"Capacity" in the graphs). 

B. Results for Horizontal QPSK 

The simulation results show that the standard deviation 
of the error is approximately 2-3dB for each stream using 
both the full-sum and max-log approximations. Moreover, the 
standard deviation of the error when using the union bound 
approximation is some ldB, which means only ldB better than 
the former approximations. These differences are well shown 
in Figure [5] We note that each graph was compensated with 
its mean error so it is symmetric around zero. It is easy to see 
that the max-log approximation based SNR and the full-sum 
approximation based SNR are very similar for stream 0. (The 
same holds for stream 1, too). Hence we may conclude that 
the max-log approximation introduces negligible performance 
degradation. 

1) Vertical QPSK: In order to expand the algorithm to 
support a single SNR we defined the SNR to be the logarithmic 
average of the 2 per stream SNRs 



SNR yertJca/ = i(SNR + SNRi), 



(35) 



for the union bound, full-sum and max-log approximations. 
Fig. [6] shows the error distribution of the vertical SNR 
calculation along with the capacity-based SNR as described 
in (0. The graph shows that there is a slight advantage 
for the new algorithm (for both the full-sum and max-log 
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Comparison of SNR Error for the 1 st Stream 




Error (dB) 

Fig. 5. Distribution of the SNR estimation error respective to stream 0. 

approximations, which similarly to the horizontal case have 
very close performance) over the capacity based SNR in terms 
of error standard deviation. As expected, the union bound SNR 
is superior to all three other calculations, having a standard 
deviation of approximately ldB. 

Error distribution of Approximated SNR for the Vertical Case 




Fig. 6. Distribution of the SNR estimation error for the vertical case 

C. 16QAM and 64QAM 

As mentioned, the SNR should not be constellation depen- 
dent. Therefore we examine the new method when 16QAM 
and 64QAM constellations are applied. The only modification 
needed is the calculation of the reference SNR which for 
QPSK was introduced in ( fT2b . In the case of 16QAM the 
error probability is given by 

SNR, 

"^0~, (36) 



p(error in stream i) = e 
and in the case of 64QAM it reads 



p(error in stream i) 



SNR, 

42 



(37) 



The values in the exponent denominators are obtained from 
the normalized constellation construction as depicted in the 
IEEE 802.16 standard 0. 



1) Preliminary verification: A straight forward extension of 
the ideas demonstrated in Section lTlI-Al is to develop a 16QAM 
based SNR estimation mechanism for 16QAM transmission 
and a 64QAM mechanism for 64QAM transmission. The 
performance of a 16QAM based mechanism for 16QAM 
transmission is given in Fig. [7] Simulation results show that 
the performance in the 16QAM case is similar to that in QPSK 
(the standard deviation is approximately 2dB) . 

Bearing in mind that SNR should be a modulation inde- 
pendent measure, we suggest to employ a QPSK based SNR 
estimation mechanism for other modulations as 16QAM and 
64QAM. Obviously, this approach leads to a more efficient 
mechanism since it implies a single mechanism for all modu- 
lations and more important, it makes use of the smaller search 
spaces respective to QPSK (the 16QAM based set, which is 
not introduced in this work, consists of 50 entries compared 
to only 13 entries in the QPSK based set, as given in d29l ). 
The performance of a QPSK based mechanism for 16QAM 
transmission is given in Fig. [7] Note that the performance 
of the QPSK based mechanism is slightly inferior to the 
16QAM based mechanism, however, this degradations seems 
negligible. 

Error distribution of Approximated SNR: QPSK vs. 1 6QAM 




Fig. 7. Distribution of the SNR estimation for QPSK and 16QAM 

2) Per Stream SNR for 16QAM: Simulation report indicated 
a standard deviation of 2.2dB for the union bound SNR 
calculation and 3dB for the approximated SNR calculation 
For the horizontal case. As for the vertical case, similarly 
to QPSK, a small degradation in means of error standard 
deviation (~ O.ldB) was achieved using the approximated 
SNR over the traditional capacity based SNR. Figure [8] shows 
the distributions of the SNR error for stream 0. 

V. Discussion and Conclusions 

In this paper, we develop methods for estimating the post 
processing SNR in ML decoded SM. Our results include a 
series of SNR estimation methods based on various approx- 
imations for the per stream error probability in ML decoded 
SM. We propose a very low complexity implementation of 
the SNR estimation method, based on the ML decoder itself 
with negligible overhead over the routine employment of 
the decoder for data detection. We show that QPSK based 
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1 6QAM : Comparison of SNR Error for the 1 st Stream 




Error (dB) 

Fig. 8. Distribution of the SNR estimation error of stream for 16QAM 
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algorithm provides plausible performance also for the higher 
modulations, so that a single SNR estimation mechanism is 
required for link adaptation. Per stream SNR estimation for 
ML decoded SM may play an important role in UL and DL 
SDMA schemes where each stream corresponds to a different 
user, hence per stream SNR estimation is required [17], |18|, 
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